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Abstract. 

We discuss a method to obtain the true event- by- event net-charge multiplicity 
distributions from a corresponding measured distribution which is subjected to detector 
effects such as finite particle counting efficiency. The approach is based on the Bayes 
method for unfolding of distributions. We are able to faithfully unfold back the 
measured distributions to match with their corresponding true distributions obtained 
for a widely varying underlying particle production mechanism, beam energy and 
collision centrality. Particularly the mean, variance, skewness, kurtosis, their products 
and ratios of net-charge distributions from the event generators are shown to be 
successfully unfolded from the measured distributions constructed to mimic a real 
experimental distribution. We demonstrate the necessity to account for detector 
effects before associating the higher moments of net-charge distributions with physical 
quantities or phenomena. The advantage of this approach being that one need not 
construct new observable to cancel out detector effects which loose their ability to be 
connected to physical quantities calculable in standard theories. 

PACS numbers: 25.75.Gz,12.38.Mh,21.65.Qr,25.75.-q,25.75.Nq 
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1. Introduction 

Higher moments of the event-by-event distribution of conserved quantities like net- 
charge, net-baryon number and net-strangeness in heavy-ion collisions have been found 
to be useful observable to characterize the system formed in the collisions pQ. Higher 
moments have been shown to be related to the correlation length pj and susceptibilities 
[3j H] of the system and hence can be used to look for signals of phase transition and 
critical point (5j EJ Ej. They have also been shown to be useful for studying the bulk 
QCD thermodynamics at high temperature [8]. Specifically, proposals have been made 
to extract the freeze-out properties of the system using higher moments of net-charge 
and net-baryon number distributions, in a way very similar to that done using the 
particle yields and ratios [HI [10J [EE] • 

Any experimental measurement is susceptible to the effects such as the finite 
acceptance, finite efficiency of counting the number of particles produced in the collisions 
and other background effects [12]. It is almost impossible to know some of these 
quantities for each event so as to correct for the effects in an event-by-event distribution. 
Hence most of the experimentally measured event-by-event distributions are presented 
without these corrections [U [131 El 03] • These corrections are carried out on an average 
level for reporting the yields of the produced particles (typically the first moment of the 
multiplicity distributions) [12J. Comparison of uncorrected experimental event-by-event 
distributions to theoretical calculations needs to be done carefully. For example, using 
the corrected mean multiplicities to explain the uncorrected measured event-by-event 
distributions could lead to different conclusions [T6| [T7]. 

Judicious construction of event-by-event observables have been proposed to cancel 
out detector effects to first order [T8| \19\ 1201 f2Tj . However, while making these 
constructs, one may sometimes loose the ability to compare them to the theoretically 
calculated quantities in order to extract meaningful physical insights. That introduce 
additional complexities which makes it difficult for a proper physical interpretation 
of the observable. As an example, the moments of the multiplicity distribution of 
conserved quantities can be shown to be proportional to correlation length (£) of the 
system. The variances (cr 2 = ((AN) 2 ); AN = N — M; M is the mean) of these 
distributions are related to £ as a 2 ~ £ 2 , skewness (S = ((AN) 3 ) /a 3 ) goes as £ 4,5 and 
kurtosis (k = [((AiV) 4 ) /cr 4 ] - 3) goes as £ 7 [2]. Their product such as kg 2 are related 
to the ratio of fourth order (x ) to second order (x^) susceptibilities [31 Hj. Where 
^,(2) _ (t^p ) ; y[ s th e volume, and AiV could be the net-baryon number or net-charge 
number. In order to cancel out the acceptance and efficiency effects to first order for 
these observables, constructs such as normalized factorial moments (defined later) can 
be made. The factorial moments of a particular order however become complicated 
function of lower order moments. Thereby making their interpretation difficult in terms 
of physical observables such as £ or x calculated in a standard theory 

Here we give a simple calculation to illustrate our point of view. Let N represents 
the produced multiplicity and n being the actually measured multiplicity in an 
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experiment. We parametrize the detector response in the experiment by a binomial 
probability distribution function given by, 

where e is the particle counting efficiency 

We further consider that the produced multiplicity follows probability distribution 
function P(N), and that for measured distribution is P(n). Then the mean of measured 
multiplicity distribution (n) can be related to the mean of the actually produced 
multiplicity distribution as, 

(n) = J nP(n)dn = J ndn J B(n | N)P(N)dN (2) 

= J P(N)dN J B{n | N)ndn = e J P(N)NdN = e(N). 

Similarly it can be shown that, 

<n 2 >= e (l- e )(A> + e 2 (A 2 >. (3) 

Now let us suppose that we can correct event-by-event particle counting efficiency, 
the variance of the resultant measured distribution can be shown to be, 

a\n/e)= l —^{N)+a 2 {N) (4) 

We find that the variance of n/e is not equal to the variance of A even though the mean 
of n/e is equal to the mean of A. Similar derivations and conclusions can be done for 
higher order moments. 

Alternatively, one can construct second order factorial moments such as 
(n(n-l)) = e 2 (A(A-l)) = (A(A - 1)) 

(n) 2 e 2 (N} 2 (iV) 2 ' 1 ' 

and the fourth order factorial moment as, 

(n(n-l)(n-2)(n-3)) = (A(A - 1)(A - 2)(A - 3)) 

(n) 4 (A) 4 ' K) 

These are found to be independent of efficiency effects. In these we also assume that e 
does not vary event-by-event. 

However a closer look at these construct will reveal that, 
(n(n - 1)) a 2 1 



(n) 2 (n) 2 (n) 



+ 1 (7) 
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While trying to remove the detector effects we have arrived at constructs which loose 
the purity of moments or become involved functions of lower order moments. Thereby 
making it difficult to directly connect to physical observables such as susceptibilities 
or their ratios which can give important insight to the bulk properties of the system 
formed in heavy-ion collisions. 

Keeping in view the importance of higher moments of multiplicity distributions to 
characterize the system formed in heavy-ion collisions, it is necessary to have a proper 
way to compare the measurements and theory calculations. At the same time ensure 
that experimental artifacts like acceptance and particle counting efficiency are removed. 
In this paper, we propose an approach based on unfolding of the measured (actually 
measured in experiments) multiplicity distribution to get back the true (actually 
produced in the collisions) distributions produced in the collisions. Such a method seems 
to work only if the detector response can be satisfactorily modeled and the statistics is 
large enough. 

The paper is organized as follows. In the next section we discuss the event 
generators used in this study. In section 3 we discuss method of unfolding. In section 
4 we present the results for the moments of the net-charge distribution as a function 
of collision centrality (defined in terms of the number of participating nucleons, N part ). 
A brief discussion of the limitations of the approach is also presented. Finally we 
summarize our study in section 5. 

2. Event generators 

In this study we have used two event generators HIJING [22] (version 1.37) and 
THERMINATOR [23j (version 2.0). They provide the possibility of different 
probability distribution for charged particle multiplicity, to study our proposal. While 
HIJING distributions are based on the physics due to QCD inspired models for 
multiple jet production, the THERMINATOR distributions are based on systems in 
thermodynamical equilibrium. The further details of the models can be found in 
Ref. [22] for HIJING and in Ref. [23] for THERMINATOR. For HIJING the events 
were generated with default settings and jet quenching on, while for THERMINATOR 
the default settings were used. We only focus on the net-charge distributions within a 
realistic acceptance of the current experiments at RHIC, that is pseudorapidity range 
between -0.5 < rj < 0.5, transverse momentum range between 0.2 < px < 2.0 GeV/c with 
full azimuthal coverage. The analysis is carried out for 19.6 GeV Au+Au collisions using 
the events from HIJING model and 200 GeV Au+Au collisions using the events from 
the THERMINATOR model as a function of collision centrality. About 5 million events 
are produced for each centrality studied in both the event generators. We have checked 
that the conclusions from each model at other energies are similar to that presented in 
this paper. Such a combination of model and beam energy is an arbitrary choice done 
to reflect a wide range of kinematics and physics of particle production. The average 
charge particle multiplicity counting efficiency is taken to be 65% following the efficiency 
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as a function of pt available for charged pions in Ref. 
3. Bayes method for unfolding of distributions 

The Bayes unfolding algorithm of RooUnfold package is used in general to remove 
the effects of measurement resolutions, systematic biases and detection efficiency to 
determine the true distributions (23]. The RooUnfoldBayes algorithm based on Bayes 
theorem uses the method described by D'Agostini in Ref. [25] . 

The procedure of Bayes unfolding can be explained by the causes C and effects 
E. In our study, causes correspond to the true multiplicity values and effects to the 
measured multiplicity values which are affected by the inefficiencies. If one observes 
n(E) events with effect E due to several independent causes (C^i = l,2,...,n c ) then 
the expected number of events assignable to each of the causes is given by: 

h(C t )=n(E)P(C t \E) (9) 

where 

p(c-\e) - (10) 

Now if we observe that the outcome of a measurement has several possible effects 
Ej(j = 1, 2, 3, He) for a given cause Cj then the expected number of events to be 
assigned to each of the causes and only due to the observed events can be calculated to 
each effect by: 

n E 

h{C i ) = Y.n{E ] )P{Q\E ] ). (11) 
i=i 

P(Ci\Ej) is the probability that different causes Cj were responsible for the observed 
effect Ej and is calculated by Bayes theorem as: 

p(am _ P(E 3 \ Cl )P (C t ) 

P{C%lE ' } ~ ^PWamct) (12) 

where -Po(Cj) are the initial probabilities. If we take into account the inefficiency then 
the best estimate of the true number of events is given by, 

i ri E 

n(Ci) = - Y."i l \ l )l'iC; /•:,) e^O (13) 

e S=i 

where is the efficiency of detecting the cause in any of the possible effects. If ej=0 
then n(Ci) is set to zero, since the experiment is not sensitive to the cause Cj. 

The above equation can be written in terms of unfolding or response matrix Mij 

as, 

n E 

n(C l ) = J2M lJ n(E J ) (14) 

3=1 
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The response matrix is constructed by repeated application of Bayes theorem and 
the regularization is achieved by stopping the iterations before reaching the "true" 
inverse. Further details of the procedure can be found in [25J. 

For the present study 5M Au+Au collision events are produced for each centrality 
bin at V /I^ = 19.6 GeV and 200 GeV using HIJING and THERMINATOR event 
generators respectively. With these events, the true distribution of net-charge (AN = 
N + — N~) is constructed on an event-by-event basis. The positive (N + ) and negative 
(N~) charge particles are selected for each event with transverse momentum range 
between 0.2 < p T < 2.0 GeV/c and pseudorapidity range between -0.5 < rj < 0.5 with 
full azimuthal coverage. 

The individual true N + and N~ are smeared with a Gaussian function with 
the mean value corresponding to the average efficiency of 65% as obtained from 
the parametrization of the pt dependent efficiency for charged pions from STAR 
experiment [12J. The width of the Gaussian distribution is taken as 10% of the mean. 
The smeared N + and N~ distributions will be called as measured distributions. The 
measured net-charge distribution is then constructed with these measured N + and N~~ 
distributions. 
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Figure 1. (Color online) Top panel: Event-by-event distribution of positive, negative 
and net-charge (denoted as "True", solid line) in Au+Au collisions for impact 
parameter b = 1.7 fm at ^/.snn = 19.6 GeV from HIJING event generator. Also 
shown are the corresponding distributions after applying acceptance and efficiency 
effects as discussed in the text (denoted as "Measured", open circles). The unfolded 
distributions are shown as red stars and denoted as "Bayes". Bottom panel: Shows 
the ratio of the unfolded to the True distributions. 



To construct the response matrix for each centrality, 2.5M events are used as 
training true distribution of net-charge and rest of the events are used as training 
measured (after smearing on an event-by-event basis) distribution. The events for 
training true and training measured are selected separately to construct the response 
matrix, in order to avoid the effect of auto-correlation. It also uses the information 
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of an event that is not measured out of true distributions and is counted towards the 
inefficiency while constructing the response matrix. 

The measured distribution of net-charge from the remaining 2.5M events is unfolded 
with response matrix obtained from the training procedure using iterative Bayes 
theorem. The number of iterations is called the regularization parameter. The present 
study uses the optimal value of 4 for the regularization. True, measured and unfolding 
are done for finer bins of each centrality and then combined to make 5% bin to eliminate 
the finite centrality bin-width effect. The moments of net-charge distributions are 
derived using cumulant method as described in [26] and are compared for true, measured 
and unfolded distributions. 

4. Results and Discussions 

Figure [H shows the true, measured and unfolded distributions for positive charge 
(panel a), negative charge (panel b) and net-charge (panel c) for most central events 
corresponding to an average impact parameter of 1.7 fm of Au+Au collisions from 
HIJING at ^/snn — 19.6 GeV on an event-by-event basis. The true distributions are 
shown as solid lines, measured distributions (subjected to particle counting efficiency) 
are shown as blue open circles and the unfolded distributions denoted as "Bayes" are 
shown as red star. For all the cases, the respective true distributions are reproduced 
from the measured distribution using the unfolding technique. The bottom panel of 
Figfj] shows the ratio of unfolded to the true distributions corresponding to the same 
distributions as shown in the respective top panels of FigHJ 

The ratio is close to unity within the statistical errors, suggesting that the unfolding 
procedure is able to get back the true distribution from a measured distribution which 
is subjected to inefficiencies in particle counting. Similar conclusions are obtained for 
distributions for Au+Au collisions at v /i^ = 200 GeV from THERMINATOR, hence 
are not shown in this paper. From now onwards we will only concentrate on the net- 
charge distributions. 

The four moments M, a, S and k of the net-charge distributions in Au+Au 
collisions at ^/snn = 19.6 GeV from the constructed true, measured and unfolded 
distributions as a function of centrality (iVp art ) are shown in FigfJJ The mean and 
standard deviation increases with A^ part , while the skewness and kurtosis decreases with 
^part- This is in accordance with the central limit theorem pp. The mean and variance 
of the measured distribution are smaller compared to those of the true, as we have 
particle counting inefficiencies for the measured case. The unfolded moments are found 
to closely follow the corresponding values of their respective true distributions. This 
can be more clearly seen from the ratio plots in FigJHJ The value of the ratio of 
unfolded to true distribution as a function of iVp art is around unity for all the four 
moments studied. Thus the unfolding method followed in this paper reproduces all the 
moments of true distribution from the measured distribution. Although not shown here, 
similar conclusions are obtained separately for the positive and negative charge particle 
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Figure 2. (Color online) Mean, standard deviation, skewness and kurtosis of net- 
charge distribution in Au+Au collisions at v / 5nn = 19-6 GeV from HIJING event 
generator. Results are shown for the True, measured and Bayes unfolded distributions 
as a function of Ar part . 



2.0 



-:i.5 



m 5 



1.0 



CO 

CC0.5 



: (a) 


I I I _ 

(N + -N) " 
AuAu 19.6 GeV : 


2.0 

m 
a 


: (b) 


HIJING 




DCX>G O O O 


m 

?1.0 
o 




exdooo © e ; 




, , T 1 T , , , 1 T , T , 1 , 


■5 

"=0.5 







100 200 300 400 



N, 



'part 



100 200 300 400 

N„ art 




100 200 300 400 



Figure 3. Ratio of the unfolded to the True net-charge distribution moments of Fig [2] 
as a function of iVp art . 
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multiplicity distributions. 

The centrality dependence of ratio of moments (cr 2 /M) and product of moments 
(Sa and ko~ 2 ) are shown in FigJU The importance of unfolding is clearly demonstrated 
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Figure 4. (Color online) Ratio (panel a) and product of moments (panel (b) and (c)) of 
net-charge distributions in Au+Au collisions at ,/snn = 19.6 GeV from HIJING event 
generator. The results are for the True, measured and Bayes unfolded distributions as 
a function of iVpart- 
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Figure 5. (Color online) Mean, standard deviation, skewness and kurtosis of net- 
charge distribution in Au+Au collisions at y^NN = 200 GeV from THERMINATOR 
event generator. Results are shown for the True, measured and Bayes unfolded 
distributions as a function of iVpart . 



by looking at the dependences of the ratio and product of moments on the iVp art . 
While for the true distribution the product of moments decreases with iV part , those for 
measured actually has an opposite trend. Indicating any physics conclusions associated 
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with variation of Sa and na 2 with N pSLTt for net-charge distributions could be highly 
misleading. However, very nice agreement between true and unfolded distributions 
are observed. They are nicely consistent even for the product of higher moments (Sa 
and na 2 ) which are very sensitive to the shape of the distributions. Suggesting that 
the unfolded distributions are well reproduced as the true distributions by using Bayes 
unfolding algorithm. 

In order to validate the applicability of unfolding algorithm for different physics 
processes, a thermal model based THERMINATOR event generator is also used. 
Figure [5] shows the centrality dependence of various moments of net-charge distribution 
in Au+Au collisions at — 200 GeV from the true, measured and unfolded 

distributions. The trends of the moments as a function of iVp art is similar to that 
seen for HIJING (Fig. |2]), although the magnitude of the moments are different. All the 
four moments of the unfolded distributions are well reproduce as the true distributions. 

Figure [6] shows the a 2 /M, Sa and na 2 as a function of iVp art of net-charge 
distributions from the true, measured and unfolded distributions. Here also, as was 
seen for the HIJING results (Fig. HJ), the ratio and products of moments from unfolded 
distributions are reproduced as true distributions up to a good extent. This suggests 
that the method proposed in this paper works equally well for parent distributions 
produced from very different particle production mechanisms as well as over a wide 
range of beam energies. 
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Figure 6. (Color online) Product of moments of net-charge distributions in Au+Au 
collisions at y/s^ = 200 GeV from THERMINATOR event generator. The results 
are for the True, measured and Bayes unfolded distributions as a function of iVp art . 



To study the effect of variation of efficiency on ratio and product of moments, the 
true distributions are smeared with a constant efficiency of 65% to obtain the measured 
distributions. Figure [7] and Fig [8] show the a 2 /M, Sa and na 2 as a function of A part 
of net-charge distributions from the true, measured and unfolded distributions with 
constant efficiency in Au+Au collisions at ^/snn = 19.6 GeV and 200 GeV from HIJING 
and THERMINATOR event generators, respectively. Panel (a) of Fig [7] and Fig |8] shows 
similar effect as for event-by-event variation of efficiency (panel (a) of Fig H] and Fig [6]) 
on the a 2 /M of the measured distributions. The effect of constant efficiency on Sa 
and ko 2 (panel (b) and (c) of Fig [JJ and Fig [8]) of measured distributions is small as 
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Figure 7. (Color online) Product of moments of net-charge distributions in Au+Au 
collisions at ^/snn = 19.6 GeV from HIJING event generator with constant efficiency 
of 65%. The results are for the True, measured and Bayes unfolded distributions as a 
function of N, 
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Figure 8. (Color online) Product of moments of net-charge distributions in Au+Au 
collisions at ^/snn = 200 GeV from THERMINATOR event generator with constant 
efficiency of 65%. The results are for the True, measured and Bayes unfolded 
distributions as a function of iVpart . 



In order to see the effect of energy dependence on our results we have carried out 
this study for net-charge distributions at midrapidity in 0-5% central Au+Au collisions 
in HIJING model for = 19.6, 27, 39, 62.4, 130 and 200 GeV. The efficiency varies 

event-by-event as per a Gaussian distribution with mean of 65% and width of 10% of 
the mean. Figure [9] shows the Mean, standard deviation, skewness and kurtosis for the 
above system as a function of beam energy. The mean and variance of the measured 
distribution are smaller compared to those of the true, as we have seen for the centrality 
dependence study (Fig. [2] and Fig. EJ). This is due the particle counting inefficiencies for 
the measured case. The unfolded moments are found to closely follow the corresponding 
values of their respective true distributions. Figure HU1 shows the cr 2 /M, Scr and na 2 
as a function of a/snn of 0-5% Au+Au collisions net-charge distributions from the 
true, measured and unfolded distributions. Here also, as was seen for the centrality 
dependence results, the ratio and products of moments from unfolded distributions are 
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Figure 9. (Color online) Mean, standard deviation, skewness and kurtosis of net- 
charge distribution in 0-5% Au+Au collisions as a function of ^/snn from HJING event 
generator. Results are shown for the True, measured and Bayes unfolded distributions 



reproduced as true distributions up to a good extent. This suggests that the method 
proposed in this paper works for parent distributions produced over a wide range of 
beam energies. 

Our study shows that it is important to correct for event-by-event detector related 
effects before proper conclusions can be obtained from higher moments studies in heavy- 
ion collisions. We have provided a method of obtaining the true distributions through 
an unfolding technique. Such a method keeps the observables same and hence has the 
advantage of being used to compare to standard theory calculations. Although this 
procedure can be easily adapted to experimentally measured distributions, it has two 
important drawbacks. Unlike the current case, where we have used an event generator 
for the study and the true distribution is available for comparison, in a real experiment 
the true distribution is unknown. Hence it is very crucial that a realistic modeling of the 
detector response and particle production is available to obtain the response matrix for 
the unfolding calculations. In most cases, the modeling of the particle production and 
the detector response is highly event generator dependent and on how realistically the 
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Figure 10. (Color online) Product of moments of net-charge distribution in 0-5% 
Au+Au collisions as a function of v /snn from HIJING event generator. The results 
are for the True, measured and Bayes unfolded distributions. 



experimental conditions are simulated. The other disadvantage is, that the procedure 
works well for high event statistics as well as high average particle multiplicity per event. 
We have seen that large uncertainties enter into the unfolded distributions if we carry 
out this study with net-protons. 

5. Summary 

In summary, we have discussed a method to obtain the event-by-event true distributions 
of net-charge from the corresponding measured distributions which are subjected to 
detector effects like finite particle counting efficiencies. The approach used is based 
on the Bayes method for unfolding of distributions. We have used event generators 
HIJING and THERMINATOR to simulate the charged particle distributions produced 
in Au+Au collisions at ^/snn = 19.6 to 200 GeV respectively. The charge particle 
counting efficiency was varied by smearing the true distributions, on event-by-event 
basis using a Gaussian function with mean 0.65 and width 0.065, to construct the 
measured distributions. We have shown that the unfolded distribution has similar mean, 
variance, skewness and kurtosis as the true distributions for all the collision centralities 
studied. The product of the moments o~ 2 /M, Scr and ko 2 which show an opposite 
trends versus iVp art for the measured distributions compared to the true distributions 
are faithfully unfolded back to agree with the true distributions. For cases where the 
efficiency of charged particle counting is constant for all events, the differences between 
the measured and true are small for the product So and na 2 compared to the ratio 
o 2 /M. The unfolding process is demonstrated to work for distributions obtained from 
widely differing physical mechanism for production of charge particle and over a wide 
range of beam energy and collision centrality. They also work for the cases where the 
charged particle counting efficiencies vary event-by-event as well as for the case where 
the efficiencies are constant. 

This method has some limitations, in terms of need for a proper modeling of the 
detector response and works well for high multiplicity and high event statistics dataset. 
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However the main advantage of this method is that we do not have to construct new 
observables which cancels out the detector effects. As the new constructs are usually 
subjected to difficulties in physical interpretation and cannot be directly compared to 
standard theoretical calculations. 
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